Draft version February 1, 2008 

Preprint typeset using LAT^X style emulatcapj v. 04/03/99 



(N 
O 
O 
(N 

X> 

IX, 

vo 
(N 



> 

oo 

(N 
O 
(N 
O 

^ ' 

6 ■ 
a: 



13 



THE OPTICAL SPECTRUM OF THE SN 1006 SUPERNOVA REMNANT REVISITED 
Parviz Ghavamian 1 , P. Frank Winkler 2,5 , John C. Raymond 3 and Knox S. Long 4 ' 5 

Submitted January 31, 2002 

ABSTRACT 

We present the deepest optical spectrum acquired to date of Balmer-dominated shocks in the NW 
rim of SN 1006. We detect the broad and narrow components of Ha, H/3 and H7 and report the first 
detection of the He I A 6678 emission line in this supernova remnant. We may have detected, at the 1.5a 
level, faint He II A4686 emission. We measure a full width half maximum of 2290 ±80 kms -1 in the 
broad component Ha line, with broad-to-narrow flux ratios of 0.84^q 01 an d 0.93j^Q '\g in Ha and H/3, 
respectively. To match these observations, our nonradiative shock models require a low degree of electron- 
proton equilibration at the shock front, T e /T p < 0.07, and a shock speed of 2890±100 kms -1 . These 
results agree well with an earlier analysis of ultraviolet lines from SN 1006. The He I/Ha and He I/He II 
flux ratios also indicate low equilibration. Furthermore, our models match the observations for mostly 
ionized (~ 90%) preshock H and mostly neutral (;> 70%) preshock He, respectively. We conclude that 
the high H ionization fraction cannot be explained by either photoionization from the reverse shock or 
relic ionization from EUV photons released in the 1006 A.D. supernova. The most plausible explanation 
appears to be photoionization from the Galactic Lyman continuum. 

Subject headings: ISM: supernova remnants: individual (SN 1006)-ISM: kinematics and dynamics, 
shock waves 



Balmer-dominated filaments (van den Bergh 1976, Long, 
Blair & van den Bergh 1988, Winkler & Long 1997). 
These filaments emit an optical spectrum dominated by 
the Balmer lines of H, produced in the same ionization 
zone as the ultraviolet lines. Balmer line profiles con- 
sist of two components: a narrow component produced 
when cold ambient neutrals are overrun by the shock and 
a broad component produced when charge exchange with 
protons generates a hot neutral population behind the 
shock (Chevalier & Raymond 1978, Chevalier, Kirshner & 
Raymond 1980). The broad component width is directly 
proportional to the postshock proton temperature, while 
the ratio of broad to narrow flux is sensitive to the electron- 
proton temperature equilibration and the preshock neu- 
tral fraction. These diagnostics are invaluable tools for 
measuring shock speeds in supernova remnants and study- 
ing the efficiency of collisionless heating processes in high 
Mach number, nonradiative shocks. 

2. SPECTROSCOPIC OBSERVATIONS 

We obtained a deep long-slit spectrum of the NW fila- 
ment in SN 1006 from the 4-m Blanco telescope at CTIO 
on 1998 June 24 (UT). The RC Spectrograph was used 
with the 527 line mm -1 grating KPGL3, the Blue Air 
Schmidt camera, and the Loral 3K CCD. This combination 
covered the spectral range 3760—7480 A at a dispersion of 
1.22 A pixel" 1 . For our slit width of 2'.'0, the instru- 
mental resolution was approximately 4.5 A (205 kms" 1 at 
Ha) along the entire spectrum. The plate scale for this 
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1. INTRODUCTION 

Among the known galactic supernova remnants (SNRs), 
SN 1006 stands out as one of the few with low extinc- 
tion (E(B - V) =0.11, Schweizer & Middleditch 1980), 
a well determined age and a reasonably well known dis- 
tance (D ~2 kpc, Kirshner, Winkler & Chevalier 1987 
(hereafter KWC87), Long et al. 1988, Laming et al. 1996, 
Winkler & Long 1997). This SNR was probably produced 
by a Type la explosion (Schaefer 1996) and appears as a 

limb brightened shell 30 in diameter at X-ray (Willingale 
et al. 1996, Winkler & Long 1997) and radio (Reynolds & 
Gilmore 1986) wavelengths. Using the Hopkins Ultraviolet 
Telescope (HUT), Raymond et al. (1995) first detected ul- 
traviolet emission lines of He II A1640, C IV AA1548, 1551, 
N V AA1239, 1243 and O VI AA1032, 1038 in SN 1006. 
Raymond et al. (1995) observed equal velocity widths for 
the emission lines detected by HUT. Combining this in- 
formation with the relative strengths of the He II, C IV, 
N V and O VI lines, Laming et al. (1996) concluded that 
the lines observed by HUT were due to collisional excita- 
tion in 2600 kms -1 nonradiative shocks with little or no 
electron- ion temperature equilibration at the shock front. 
In this case the high shock speed and low density (n <i 
1 cm -3 ) result in negligible radiative cooling downstream 
(hence the term 'nonradiative shock'). All of the line emis- 
sion is produced in a narrow ionization zone close to the 
shock front . 

In the optical, SN 1006 appears as a faint network of 
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configuration was 0'.'51 pixel 1 and the slit length was 5 . 
We rotated the slit to PA = 63.9° and placed its center at 

coordinates 15 h 02 m 13!5, -41°45'22" (2000), as shown in 
Fig. 1. We chose this position to cut obliquely across the 
NW filament in order to capture as much of the bright re- 
gion as possible while giving enough background sky, away 
from the filament, for a good sky subtraction. It is virtu- 
ally identical to the position used by KWC87. We obtained 
4 exposures at this location, with a total exposure time of 
8400 s, which we combined into a single two-dimensional 
spectrum following initial reduction. 

The data were reduced using standard IRAF 6 proce- 
dures of bias subtraction, flat fielding and illumination 
correction. Spectra of a HeNeAr lamp taken at the begin- 
ning and end of the series of object frames were used for 
wavelength calibration. We subtracted the night sky con- 
tribution from the two-dimensional spectrum using emis- 
sion from the two ends of the slit. The Ha line profiles 
(both broad and narrow components) were nearly constant 
along the observed filament; therefore, we integrated the 
emission along the slit to obtain a single high signal-to- 

noise spectrum from a 51 section of the filament (Figs. 2 
and 3). Exposures of several spectrophotometric standard 
stars from Hamuy et al. (1994) were used for flux cali- 
bration. We estimate an absolute photometric accuracy of 
20% for our quoted emission line fluxes. 

3. DETECTED EMISSION LINES 

Among the most interesting new features in our optical 
data are the detection of broad H/3 and He I A 6678. To our 
knowledge, this is the first detection of the A6678 line in 
a pure nonradiative shock. The neutral He atoms passing 
downstream are unaffected by the electromagnetic turbu- 
lence at the shock front. Therefore, neutral He atoms re- 
main cool throughout the shock and produce narrow emis- 
sion lines unresolved in our data. On the other hand since 
Hc + is an ion, it is heated at the shock front by the same 
collisionless processes which heat the electrons and ions. 
Therefore, we expect the He II A4686 line to be broad and 
more difficult to detect than the He I A6678 line. There 
is a hint of the He II A4686 line at a low (1.5c) statistical 
significance. The surface brightnesses of the detected emis- 
sion lines from NW SN 1006 appear in Table 1. We have 
used E(B - V) = 0.11 (Schweizer & Middleditch 1980) 
to obtained dereddened line intensities from our measured 
spectra. 

4. SPECTRAL LINE FITS 

We fit the line profiles using a standard nonlinear \ 2 
reduction routine. We assumed Gaussian profile shapes 
for all lines and fitted the Ha, He I, H/3, H7 and He II 
lines separately. Each Balmer line profile was fit with two 
Gaussians, one each for the narrow and broad components. 
The He lines were fit with single Gaussians. To reduce the 
required degrees of freedom, we set the widths of the nar- 
row component Balmer lines and He I line equal to the 
instrumental resolution during the fitting process. 

Due to the faintness of the A4686 emission, we cannot 
estimate the width of this line from the fits themselves. 
Therefore, we set the velocity width of the He II emission 

6 IRAF is distributed by the National Optical Astronomy Observatories 
with the National Science Foundation 



line width equal to that of the Ha line. This assumption is 
supported by the HUT analysis of Raymond et al. (1995), 
who found that to within the errors, the velocity widths 
of the He II A 1640 and broad Ly /3 lines are equal (i.e., 
electron-ion and ion-ion temperature equilibration is inef- 
ficient in SN 1006). The best fit for the He II A4686 line 
profile appears in Fig. 4. We note that due to the large 
width and faintness of the He II line, the flux predicted 
by the line fit is quite sensitive to the baseline level: the 
lower the baseline, the higher the predicted He II flux. 
Both the statistical uncertainty and the systematic base- 
line uncertainty contribute to the large error bars on the 
flux. For the purpose of our subsequent analysis, the most 
important quantity here is the upper limit on the He II 
line flux. 

Our best fit for the Ha broad component width is 
2290±80 kms -1 , agreeing very well with the value of 
2310±210 km s" 1 measured by Smith et al. (1991). The 
line center shift between the broad and narrow compo- 
nent lines is 29±18 kms -1 , less than the uncertainty in 
broad component width, indicating a nearly edge-on view- 
ing angle (| 90° — O | < 2°). The Ha broad-to- narrow ratio 
is 0.84+q'q-] 3 , lying slightly above the extreme upper limit 
quoted by Smith et al. (1991), 0.73±0.06. The H/3 broad- 
to-narrow ratio is somewhat larger, 0.93^16- The H/3 
broad-to-narrow ratio is larger than that of Ha due to 
the lower conversion efficiency of Ly 7 photons into H/? in 
the narrow component (c.f. Ghavamian et al. 2001). The 
broad-to-narrow ratio of the H7 line (0.63^0.25) is poorly 
constrained due to the faintness of the broad component. 
The observed Ha surface brightness of the Balmer filament 
in SN 1006 is (3.2±0.7) x 10~ 16 ergs cur 2 s" 1 arcsec" 2 (in- 
cluding the 20% photometric uncertainty), in reasonable 
agreement with the value quoted by KWC87. 

5. ANALYSIS AND MODELS 

To model the Balmer line profiles, we used the 
one-dimensional, plane-parallel shock code described in 
Ghavamian et al. (2001) to compute a grid of numerical 
models. Briefly, the code predicts the density and temper- 
ature of electrons, protons and hot neutrals behind a non- 
radiative shock. Broad to narrow flux ratios are computed 
by a Monte Carlo simulation which follows collisionally ex- 
cited Ly (3 and Ly 7 photons until they are absorbed or 
escape from the grid. Collisional excitation and ionization 
of H by both electrons and protons is included in the mod- 
els. Free parameters of the model are the preshock neutral 
fraction f H o, shock speed, vs, and fractional equilibration 
immediately behind the shock. The latter quantity de- 
notes the degree of electron-ion equilibration by plasma 
turbulence at the shock front, where f eq = represents no 
(pure Coulomb) equilibration, while f eq = 1 represents full 
equilibration. The electron and proton temperatures are 
then defined by the relation 
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where fj, is the mean molecular weight ( = 0.6 for Galac- 
tic abundances, Allen 1973) and ^^m p w| is the mean 
postshock temperature. For an input f eq , the subsequent 
evolution of T e and T p by Coulomb collisions is followed as 
a function of position behind the shock. We have updated 
the code to include the calculation of both the He ioniza- 
tion balance and the contribution of He ionization to the 
postshock electron density. To increase the signal-to-noise 
of our simulations, we have also increased the number of 
excitations followed by the Monte Carlo simulation from 
10,000 used in the earlier paper (Ghavamian et al. 2001) 
to 50,000. 

Using the same method outlined in Ghavamian 
et al. (2001), we computed a grid of models in equilibra- 
tion and shock velocity, combining each f eq with the vs 
required to match the observed FWHM of the broad com- 
ponent. We sampled the subset f eg = (0, 0.03, 0.05, 0.1, 
0.2, 0.4, 0.6, 0.8, 1), corresponding to i; 5 = (2865, 2880, 
2895, 2920, 2965, 3100, 3235, 3395, 3580) kms" 1 , respec- 
tively. The Ha and H/3 broad-to-narrow ratios predicted 
by our models are shown in Fig. 5. It is immediately obvi- 
ous that multiple combinations of f eq and h / I n are consis- 
tent with the observed broad-to-narrow ratios. This makes 
the unambiguous, simultaneous determination of fu° and 
f eq from Fig. 5 impossible. In an earlier calculation, Lam- 
ing (2000) obtained similar shapes for the broad-to-narrow 
ratio curves, the main difference being that our curves flat- 
ten out and then decline slightly at high equilibrations. 
The reason for the disagreement is that unlike our models, 
the Laming (2000) calculation maintains a constant shock 
speed while varying f eq . The slight downturn in our com- 
puted narrow ratio curves is caused by the decline of the 
H + — H° charge exchange rate at high equilibrations (i.e., 
high shock speeds). As a result, the hot neutral density 
peaks at progressively smaller values behind the shock as 
vs increases beyond 3000 kms" 1 . 

Fortunately, we can use the He I/He II and He I/Ha flux 
ratios to break the degeneracy in Fig. 5. Since every He 
atom in the postshock ionization zone will spend some time 
as He+ (the prcshock He ++ fraction typically amounts to 
only a few percent), the predicted He II A4686 flux is in- 
sensitive to the preshock He ionization fraction. Further- 
more, as demonstrated by Hartigan (1999) and Laming 
et al. (1996), the He II A4686 flux is also insensitive to 
the degree of electron-ion temperature equilibration be- 
hind the shock (a result that is confirmed by our models). 
Therefore, we can use the He I/He II ratio to constrain the 
prcshock He fraction. Once this quantity is known, we 
can then use the He I/Ha ratio to constrain the preshock 
H° fraction. This will enable us to determine the appropri- 
ate choice of broad-to-narrow ratio from Fig. 5. From our 
shock models we estimate that regardless of the assumed 
equilibration, 90% of the He II emission is produced within 
" 

a zone <1 thick (assuming a total prcshock density of 1 
cm" 3 (Winkler & Long 1997), shock speed 3000 kms" 1 
and a distance of 2 kpc). This is smaller than both the 
pixel scale (0'.'51 pixel" 1 ) and the slit width (2'.'0). There- 
fore, we believe that all of the He II emission layer was 
contained within the slit. 

In most interstellar shocks where He I emission is ob- 
served (mainly radiative shocks), the A5876 line is brighter 
than the A6678 line. The A6678 emission line arises from 



the singlet transition 1 P (ls2p) — 1 D (ls3d) and is most ef- 
ficiently excited from the ground state (i.e., by collisional 
excitation). The A5876 line, however, arises from a triplet 
transition and is most efficiently excited by cascades from 
high lying energy levels (i.e., from recombination). There- 
fore, since recombination is negligible in SN 1006 shocks, 
we expect the former line to be significantly brighter than 
the latter (indeed, A5876 emission is not detected in our 
optical spectrum). In the following analysis we limit our 
calculations to the A6678 line intensities. 

To estimate the He I/Ha ratio for SN 1006, we com- 
puted the total Ha and He I A6678 fluxes through the 
shock front using the combinations of f eq and v$ listed 
above. Calculation of the He I A6678 flux in SN 1006 is 
relatively straightforward because proton excitation and 
ionization of He is negligible (Laming et al. 1996). In addi- 
tion, the only significant contribution to the He I flux other 
than *S (Is 2 ) — X D (ls3d) excitation are cascades from the 
4f level via i D (ls3d) - Y F (ls4f). We have included both 
direct excitation and cascade contributions to 3d in our 
He I A6678 flux calculation. We obtained the electron ex- 
citation cross sections from the convergent close coupling 
(CCC) results of Fursa & Bray (1995). At high energies, 
we approximated the cross sections with a Born-like de- 
pendence, scaled to match the CCC values at 0.9 keV. 
Our He ionization balance calculation uses the electron 
ionization rates of Younger (1981) for He and He + . We 
also utilized the CCC cross sections of Fursa and Bray 
(1995) to compute the collisional excitation rate of He II 
A4686 electrons (proton excitation is negligible, see Lam- 
ing et al. 1996). Given the low density, the order of mag- 
nitude lower abundance of He relative to H, and the large 
line width of He II, we assumed Case A conditions in the 
He II line, i.e., the radiative excitation of He II A4686 
by He II A243 (Ly 7) absorption is negligible behind the 
shock. 

The model F(6678)/F(4686) ratios are shown in Fig. 6 
for a range of preshock He fractions and temperature 
equilibrations. The dashed lines indicate error bars on 
the observed line ratio after dereddening. Since the ra- 
tio of the He I excitation rate to the He I ionization rate 
depends steeply on electron temperature for T e <J 10 5 K, 
the He I/He II ratio declines sharply at low values of f eq 
(<0.01). It is clear from this figure alone that matching 
the predicted He I/He II ratio with the observed value 
requires a preshock He neutral fraction > 0.5 and low 
electron-ion equilibration. 

The predicted He I/Ha ratios are displayed in Fig. 7 for 
a range of preshock H neutral fractions. For simplicity, 
we assumed /jj e ° = 1 in these models (varying this quan- 
tity merely moves the curves up or down). To facilitate 
comparison with the broad-to-narrow ratios, we used the 
same set of neutral fractions and combinations (f eq , vs) 
employed by the models shown in Fig. 5. 

(1) Case A in both the broad and narrow components, 
and (2) Case A in the broad component, Case B in the nar- 
row component. Since the inferred Ly (3 and Ly 7 optical 
depths behind Balmer-dominated shocks are often <J 1 in 
the narrow component line (Chevalier, Kirshner & Ray- 
mond 1980, Ghavamian 1999, Ghavamian et al. 2001), the 
Case A and Case B line ratio curves mark the largest and 
smallest values of He I/Ha attainable for each combination 
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of preshock H and He ionization fraction. 

Comparing Figs. 6 and 7, we find that (f ) the observed 
He I/He II and He I/Ha ratios are simultaneously matched 
for f H o ~0.1 and f He o > 0.7 in the preshock gas, and 
(2) only the lowest equilibration models (f eq < 0.05) can 
match the observations. Returning to Fig. 5 and choosing 
the lb/ I n curves for fjjo =0.1, we find a match between 
the observed and modeled Ha and H/3 broad-to-narrow ra- 
tios for / eg < 0.1, or T e /T p < 0.07. This corresponds to a 
shock speed of 2890±100 kms -1 (including the measure- 
ment uncertainty of the broad component). This result 
is in good agreement with the shock speed estimated by 
Laming et al. (1996) from ultraviolet line analysis of HUT 
data: 2600±300 kins" 1 . It also agrees with the low equi- 
libration inferred by Laming et al. (1996): T e /Ti <0.05. 

Given the intrinsic faintness of the He II emission, it 
is useful to compare the measured He II A4686 surface 
brightness with the value inferred from observation of the 
He II A1640 line by HUT (Raymond et al. 1995). Assum- 
ing Case A conditions in the He II lines (reasonable due 
to the low density of He), the predicted He II A4686 flux 
in the RC Spectrograph slit is 
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where 74686 is the number of A4686 photons per in- 
cident He + ion and 71640 is the corresponding quan- 
tity for the A1640 line. We have used our Ha image 
(the star-subtracted version of Fig. 1) to scale from the 
He II A4686 flux expected in the HUT aperture to the 
value expected in the RC Spectograph slit. This scal- 
ing is valid if the H neutral fractions are roughly equal 
in the two filaments (reasonable given their close prox- 
imity and comparable surface brightnesses). From the 
Ha image we estimate Fncspec(F[a) / 'FnuT(Ha) «0.13. 
Next, from the calculations of Hartigan (2000) we esti- 
mate 74686/71640 ~ 0.1 for 10 5 K < T e < 10 7 K. Inserting 
these values into Eqn. || and using Fjf;yr( 1640) from Ray- 
mond et al. (1995, who also assumed E(B — V) =0.11), 
we find F(4686) = (1.1±0.2) x 10~ 15 ergs cm~ 2 s _1 . This 
agrees well with the A4686 flux we derive from our data, 
(9.4±6.4)xl0^ 16 ergs cur 2 s" 1 . 

6. DISCUSSION 

In their detailed imaging study of SN 1006, Winkler 
& Long (1997) found that X-ray emission in the NW 

peaks approximately 12 behind the optical filaments. Us- 
ing nonequilibrium ionization models convolved with the 
ROSAT HRI response function, they found that a total 
preshock density ~ 1 cm~ 3 was required to match the spa- 
tial distribution of the X-ray emission. Combined with 
our estimated H preshock neutral fraction, this implies a 
preshock neutral density n^o ^0.1 cm -3 with tl v p ~ 0.5 
in the narrow component. The inferred value of tl v p will 
be even lower if Winkler & Long (1997) overestimated the 
total preshock density. The dereddened narrow Balmer 
decrement is 3.35±0.14, significantly lower than the value 
of 4.3-6 observed in Tycho's SNR (Ghavamian 2001). The 
absorption of narrow component Lyman photons within 
the shock structure scales linearly with fn , as long as fn° 
is small (<; 0.5). This enhances the Ha emission more than 
H/3, steepening the narrow component Balmer decrement. 



Therefore, the modest decrement we observe in SN 1006 
is consistent with our low derived Lyman optical depth. 

6.1. The Preshock Ionization Fractions 

The detection of the He I A6678 line indicates that a 
significant fraction of the preshock He is neutral. The 
collisional excitation of He II A304 behind the shock cer- 
tainly produces some photoionization and heating of the 
ambient gas (see, for example, Tycho's SNR, Ghavamian 
et al. 2000). Since the photoionization rate of He by A304 
photons is nearly an order of magnitude higher than the 
photoionization of H° (Ghavamian et al. 2000), some addi- 
tional mechanism is required which can ionize ambient H° 
more effectively than He . There are several possibilities 
which we discuss in turn below. 

6.1.1. Photoionization From the Reverse Shock? 

One process which may generate high H ionization frac- 
tions in the ambient medium while leaving most He neu- 
tral is photoionization by O II— O VIII emission from the 
outermost layer of shocked ejecta (e.g., Hamilton & Fesen 
1988, hereafter HF88). The strongest of these lines have 
energies lying between the H and He ionization thresh- 
olds. Using the calculations of HF88 and the He II A304 
photoionization models of Ghavamian et al. (2000), we can 
estimate the ratio of the H photoionization rate to the He 
photoionization rate in the interstellar medium due to ra- 
diation from the blast wave and shocked ejecta. Since each 
O atom entering the reverse shock will pass through all 
ionization stages between ++ and +7 , the ratio of pho- 
toionizations ajj /He (R) a t a distance R from the remnant 
center is 



&H/He(R) 

where 



l(K)o- H (K) +7/7(304)^(304) 
7(A l )<r He (A l ) + 7/7(304)^(304) 



(4) 
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Here, ^(\) is the number of ionizing photons per ion 
of wavelength A; (the Bethe parameter as tabulated by 
HF88 for the reverse shock), 7(304) is the Bethe param- 
eter for He II A304 photons from the blast wave (w 1.1 
for vs >2000 kms -1 , Ghavamian et al. 2000), a H is the 
H° photoionization cross section from Osterbrock (1989) 
and a He is the corresponding cross section for He (Reil- 
man & Manson 1979). The current O-layer and blast wave 
radii are R e j (O) and Rb w , respectively, while N e j (O) is the 
number of shocked O atoms in the ejecta and Nb w (He) is 
the number of He atoms swept up by the blast wave. In 
our calculation we have assumed that the ionization state 
of the ISM is negligibly affected by recombination. HF88 
found that the relative masses of Fe, Si and O predicted 
by carbon deflagration model CDTG7 (Woosley & Weaver 
1987) matched the observed Fe II AA2586, 2599 absorption 
profile in SN 1006 best. Therefore, using the shocked O 
ejecta mass of 0.41 Mq and blast wave mass of 4.3Mq 
from model CDTG7, R e j(0) — 6 pc and Rb w =8.3 pc 
from HF88 and neglecting radiative transfer effects on the 
photoionizing radiation, we find an/He ~0.1 for a parcel 
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of interstellar gas just outside the blast wave radius. This 
result suggests that reverse shock photoionization cannot 
produce our derived low H neutral fraction. The smart- 
ness of an/He is mainly due to that fact that (1) there are 
nearly an order of magnitude more He atoms in the blast 
wave than there are O atoms in the reverse shock, and (2) 
since R e j (O) < Rb w , a given parcel of interstellar gas sees 
a larger emitting area from the blast wave than from the 
reverse shock. 

A more detailed, time dependent calculation would be 
desirable. However, we note that the reverse shock has 
already propagated through the entire O-layer, while the 
blast wave continues to accumulate interstellar material 
(and hence He atoms). Therefore, the ratio M e j(0)/Mb w 
is decreasing as the supernova remnant evolves and the 
blast wave is becoming the increasingly dominant pho- 
toionization source for the ISM around SN 1006. Finally, 
we note that even if we had included photoionization from 
the Si-rich and Fe-rich layers in Eqn. ^, our conclusions 
would remain unchanged since most of the EUV/X-ray 
lines from Si and Fe will ionize He much more effectively 
than H. 

6.1.2. Relic Ionization from the Supernova Explosion? 

An interesting source of photoionization is the EUV 
flash produced when photons escape from the detonating 
white dwarf (Type la SN). Using the EUV burst models 
presented by Blinnikov & Sorokina (2000) and Sorokina & 
Blinnikov (2001, in preparation), we estimate ~ 10 55 — 10 56 
ionizing photons (E > 13.6 eV) are released over the inte- 
grated Type la SN light curve. The precise number can 
vary by a factor ~5, depending on which detonation model 
is used in the light curve calculation (see Blinnikov & 
Sorokina (2000) for model descriptions). Due to the steep- 
ness of the EUV spectrum, the H° ionization rate is much 
larger than the He ionization rate, a desirable feature. 
However, assuming an average ISM density of 0.06—1.0 
cm~ 3 around SN 1006 and distance of 2 kpc, there were 
;> 5x 10 57 H atoms contained within the current blast wave 
radius at the time of the SN explosion. Even if these H 
atoms were already partially ionized beforehand, the total 
ionizing flux from a Type la explosion still falls an or- 
der of magnitude short of the value needed to effectively 
ionize the gas at the current blast wave radius. The dis- 
parity is worsened if we include radiative transfer loss of 
the released photons as they (1) propagate through the 
ISM, and (2) possibly encounter a circumstellar envelope 
left by a nondegenerate donor star prior to the explosion 
(Cumming et al. 1996). Therefore, regardless of whether 
the EUV emission ionizes H more effectively than He, there 
are not enough photons produced in the explosion to ex- 
plain our inferred neutral fractions. This result also agrees 
with the SN flash ionization calculations of Cheng & Bruh- 
weiler (1996). 

6.1.3. Photoionization from Galactic Background 
Radiation? 

If neither the supernova event nor the postshock emis- 
sion from the supernova remnant can account for our de- 
duced H° and He fractions, we must consider nonlocal 
sources of ionizing radiation. Given the lack of late O and 
early B stars in the vicinity of SN 1006, the most promis- 
ing source is the Galactic ionizing background. Reynolds 



(1984) observed the Ha recombination line toward a sam- 
ple of high Galactic latitude pulsars (| b \ > 5) at distances 
of 2—3 kpc. He found that Lyman continuum radiation 
from O stars in the Galactic plane was the most likely 
source of ionization in the diffuse ISM. Since the Lyman 
continuum is produced by recombination, it declines nearly 
exponentially beyond each ionization edge (Brussard & 
Van De Hulst 1962). Therefore, incident Lyman contin- 
uum photons will ionize H far more effectively than He, a 
desirable property. 

The remaining question is whether the Galactic Lyman 
continuum flux is high enough at the position of SN 1006 
(b = +14.6°) to produce the necessary amount of H ioniza- 
tion. Assuming a density of 1 cm~ 3 around SN 1006, the 
recombination rate for 5000 K^T^10,000 K is ~5xl0~ 13 
s _1 (Osterbrock 1989). Therefore, obtaining an H neutral 
fraction ~0.1 requires a photoionization rate ~5xl0~ 
s -1 . Dove, Shull & Ferrara (2000) estimated a Lyman con- 
tinuum flux ~ 5 x 10 7 photons cm~ 2 s _1 from O stars in the 
Galactic disk. In their time dependent radiative transfer 
calculation for the ionizing continuum, they estimated that 
roughly 3%— 6% of the Lyman continuum photons escape 
through superbubbles and Galactic chimneys into the halo 
of the Milky Way, producing the "Reynolds Layer" of dif- 
fuse Ha emission (Reynolds 1984). If the Lyman photons 
propagating into the halo suffer no further attenuation, the 
corresponding H ionization rate in the vicinity of SN 1006 
is approximately Fl v o- h (912) ~1— 2xl0~ n photons s^ 1 , 
more than enough to produce the necessary H ionization 
around SN 1006. 

If the Galactic Lyman continuum background pre- 
ionizes most of the ambient H on the NW side of SN 
1006, the question arises of why the ambient H on the 
eastern half of Tycho's SNR (also the remnant of a Type 
la explosion) was found to be mostly neutral (Ghavamian 
et al. 2000). First, the eastern half of Tycho's SNR is 
likely interacting with the outer edge of a warm H I cloud 
(Reynoso et al. 1999, Ghavamian et al. 2000). Second, 
the presence of large structures in the Galactic disk (H I 
clouds, superbubbles, molecular cloud complexes, etc.) 
produces large variations in the Lyman continuum flux 
from position to position and time to time (Dove, Shull & 
Ferrara 2000). Since Tycho's SNR is close to the Galac- 
tic plane (b =+1.4°), it is possible that this remnant is 
shielded from the Lyman continuum emission by the many 
dense clouds in the Galactic disk. SN 1006, on the other 
hand, lies nearly 500 pc above the Galactic plane. There- 
fore, neutral H near SN 1006 may "see" a larger time aver- 
aged Lyman continuum than the neutral H around Tycho's 
SNR. 

6.2. Comparison With Ionization Estimates of the Local 

Bubble 

An interesting comparison can be drawn between our de- 
rived H and He neutral fractions and values estimated for 
the local (<; 200 pc) ISM. For example, using EUVE spec- 
tra of 29 DA white dwarfs, Wolff et al. (1999) estimated 
fn e + ~ 0.4 in the Local Bubble, with no evidence of a gra- 
dient or other general pattern. They interpreted this result 
as evidence for a lack of He-ionizing sources in any partic- 
ular location. On the other hand, Wolff et al. (1999) found 
a positive gradient in the H ionization fraction toward the 
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Canis Major cavity, implicating the B-stars Adara and 
Mirzam as the photoionizing sources. In some respects 
the ionization properties of the Local Bubble and the ISM 
near SN 1006 are similar: the H ionization is strongly in- 
fluenced by EUV radiation, while He remains largely un- 
affected. The main difference between the two cases, of 
course, is the lack of a local H-ionizing source near SN 
1006. Our derived He ionization fraction (< 0.3) for SN 
1006 is lower than the average value of 0.4 measured by 
Wolff et al. (1999), but there are individual cases from 
their WD sample with He ionization fractions similar to 
ours. 

6.3. Comparison With Earlier Results 

Our optical analysis lends strong support to the picture 
developed by earlier X-ray, UV and optical studies of SN 
1006. Multi-epoch optical imagery of SN 1006 has yielded 

a proper motion of 0.3 ±0.04 yr (Long, Blair & van 
den Bergh 1988) for the Balmer filaments in the NW. Com- 
bined with our estimated shock speed this gives a distance 
of 2.0^0 3 kpc, agreeing well with the value 1.8±0.3 ob- 
tained by Laming et al. (1996). The SN 1006 blast wave 
is the fastest shock we have thus far modeled with our 
Balmer-dominated shock code. In modeling Balmer spec- 
tra from the Cygnus Loop (vs — 300 kms -1 ), RCW 86 
(v s -600 kms" 1 ) and Tycho's SNR (v s -2000 kms" 1 ), 
we noted that the equilibration f eq declined from near 
1 for the slowest shocks to <^0.2 for the fastest shocks 
(Ghavamian 1999, Ghavamian et al. 2001). The shocks in 
SN 1006 are closer to 3000 kms -1 and are matched by an 
even lower f eq , further evidence of an inverse correlation 
between shock speed and equilibration. The case for low 
equilibration at high Mach numbers has recently received 
further support from X-ray observations of SN 1987A. In 
their analysis of high resolution HETG data from Chan- 
dra, Michael et al. (2001) found that a high blast wave 
speed (3500 kms" 1 ) along with a low electron tempera- 
ture (f e q < 0.1) were required to match the X-ray profiles 
of the SN 1987A spectrum. 

6.4. Implications for Collisionless Shock Physics 

There arc few numerical simulations available of high 
Mach number, collisionless shocks. Cargill & Papadapou- 
los (1988) modeled the electron-ion heating in a shock 
moving perpendicular to the ambient magnetic field, with 
M A =50 [M A = v s /{c 2 s + v A fl 2 is the Alfvenic Mach 
number, where c s is the sound speed and va is the Alfven 
speed of the preshock gas) . In their hybrid simulation the 
electrons were treated as a resistive fluid, while the ions 
were treated as individual particles. Cargill & Papadapou- 
los (1988) found that two-stream plasma instabilities pro- 
duced by ions reflected upstream resulted in T e /Ti — 0.25 
at the shock front. Shimada & Hoshino (2000) tried an 
approach where both electrons and ions were treated as 
particles. At the highest Mach numbers covered by their 
simulations, they predicted that the equilibration would 
increase with Mach number. In their fastest simulated 
shock (Ma = 20), they found T e /Ti -0.2. Their curve 
of T e /Ti vs. Ma flattens considerably as M A ^20; ex- 
trapolating this curve to Ma J> 200 leads to significantly 
larger equilibration than implied by the ultraviolet (Lam- 
ing et al. 1996) and optical (this paper) observations of SN 



1006. 

During the last few years, a considerable effort has been 
made in modeling collisionless heating processes in solar 
wind shocks (M A < 20). Hull & Scudder (2000) combined 
an electron fluid treatment with constraints from ISEE 1 
data of Earth's bow shock to model the electron-ion equi- 
libration process. For the first time, they were able to 
match the observed electron temperature jumps in fast 
shocks propagating both parallel and perpendicular to the 
magnetic field. These models emphasize the role of coher- 
ent (DC) electric forces in the collisionless heating, rather 
than wave-particle interactions. The advantage of the Hull 
& Scudder (2000) treatment is that it draws upon in situ 
observations of real collisionless shocks rather than com- 
puter simulations. Both these models and the observations 
(Schwartz et al. 1988) predict T e jT l - 0.25 for shocks in 
the solar wind. Although we are unable to perform similar 
detailed observations of supernova shocks, extending the 
models of Hull & Scudder (2000) to much higher Mach 
numbers may be an interesting direction for future work 
on collisionless shocks. 

7. SUMMARY 

The optical spectrum presented here is the deepest yet 
obtained of the nonradiativc filament along the northwest 
rim of SN 1006. Lines of Ha, H/3, H 7 , RS, He I A6678 
and possibly He II A4686 have been detected— more lines 
than in the spectrum of any other purely nonradiative SNR 
filament. We have detected both broad and narrow com- 
ponent emission lines of Ha, H/3 and H7. The broad and 
narrow components result, respectively, from fast neutrals 
produced by charge exchange and slow neutrals that are 
collisionally excited. 

The broad component widths, the ratios of broad to 
narrow flux in Ha and H/3, and the Balmer line strengths 
relative to He I and He II enable us to determine the im- 
portant parameters for the shock. We find a shock speed 
of 2890±100 kms -1 and a very low degree of electron-ion 
equilibration, T e /T p < 0.07, both consistent with earlier 
measurements. In addition, our models can match the ob- 
servations only if the preshock H is mostly ionized (— 90%) 
and preshock He is mostly neutral (<;70%). In seeking 
an explanation for this ionization structure, we have con- 
cluded that neither radiation from the blast wave nor ra- 
diation from the reverse shock can effectively ionize the 
ambient H surrounding SN 1006. We also find that the 
EUV flash from the supernova explosion is incapable of 
producing the necessary H ionization. The most plausi- 
ble explanation for the high H ionization fraction seems to 
be photoionization from the Galactic Lyman continuum. 
Finally, comparing our derived electron-ion equilibration 
with values derived for other nonradiative SNRs reinforces 
the earlier conclusion that the fastest shocks have the least 
equilibration. A thorough understanding of this trend re- 
mains as a challenge for the study of collisionless shock 
physics. 
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Fig. 1. — A narrow band Ho image of NW SN 1006, acquired from the CTIO 0.9-m telescope in 1998 June, shortly before the spectroscopy 
reported here. The location and PA of the 2 X 5 RC Spectrograph slit is marked. North is at the top, East is to the left. 
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Fig. 2. — Full one-dimensional spectrum of the NW SN 1006 Balmcr filament. Narrow emission lines of Ha, H/3, H7 and H<5 are detected, 
along with broad emission lines in Ha, H/3, and possibly H7. The sharp feature near the middle of the spectrum is an artifact of the sky 
subtraction. 
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Fig. 3. — Closeup view of the NW SN 1006 Balmcr spectrum. No smoothing has been applied. Among the newly detected lines is He I 
A6678. There may also be a weak detection of the He II A4686 line. 



11 



T 




2 I I i i i i I i i i i I i i i i I 

4600 4650 4700 4750 

Wavelength (A) 

Fig. 4. — The He II A4686 line profile. The fit has been performed assuming the He II line width equals that of Ho. 
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Fig. 5. — The calculated Ha and H/3 broad-to-narrow ratios vs. the electron-ion equilibration fraction f eq at the shock front. The dashed 
horizontal lines mark the upper and lower limits of the observed broad-to-narrow ratios. Curves of constant preshock neutral fraction are 
shown for cases where the preshock H is mostly neutral, half neutral and mostly ionized. 
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Fig. 6. — The observed vs. predicted He I/Hc II line ratio, shown for a range of prcshock He neutral fractions. The observed range of the 
He I/He II ratio is indicated by the horizontal dashed lines. The model ratios have been computed assuming that neutral H in the narrow 
component is optically thin to Ly f) photons behind the shock (Case A). 
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Fig. 7. — The dependence of the He I/Ha ratio on f eq , shown for a sampling of preshock H neutral fractions. The preshock He is assumed 
to be fully neutral. The dashed lines indicate the observed range of the He I/Ha flux ratio. The ratios in the upper figure are computed 
assuming that neutral H in the narrow component is optically thick to Ly f3 behind the shock (Case B), while the ratios in the lower figure 
assume that neutral H in narrow component is optically thin to Ly fi behind the shock (Case A). 
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Table 1 

Observed Surface Brightnesses For NW SN 1006 a 



Line 


F(A) o6s 


I(A) b 


Broad to Narrow Ratio 


Velocity Width (kms^ 1 ) 


Ha 


355.3±34.8 


318.1±31.2 


84 + 03 


2290 ±80 


H/3 


100 c 


100 d 


93 + 18 


2290 (fixed) 


H 7 
H<5 

He I A6678 
He II A4686 


36.4±7.4 

>16.8 
2.0±0.4 
6.6±4.5 


38.1±7.7 

>18.3 
1.9±0.44 
7.0±4.8 


U.DO_ 25 
_f 


2290 (fixed) 
215 (fixed) 
215 (fixed) 
2290 (fixed) 



a Balmer line surface brightnesses are quoted for the combined broad and narrow component emission from a 2 x51 
section of filament. The error bars are la. 

b Dcreddened intensities assume E(B-V) = 0.11±0.02 (Schweizer & Middleditch 1980). 

c Observed total H/3 surface brightness of 9.0xl0~ 17 ergs cm~ 2 s _1 arcsec~ 2 . 

d Dereddened total H/3 surface brightness of 1.3xl0 -16 ergs cm -2 s _1 arcsec~ 2 . 

c Narrow component line widths are fixed at 205 km s -1 (the instrumental resolution). 

f The He I line is intrinsically narrow, while the He II line width is intrinsically broad. 



